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I. INTRODUCTION 


Acoustic methods have for a long time proved useful for 
deterrining the position of objects submerged in the ocean 
and continue to be useful today, especially for monitoring 
the realistic testing of advanced underwater weapons (such 
as torpedces), in the real ocean envirchment. 


One method sormmonly used requires installation of a 


‘sound source on the underwater vehicle to be tracked. The 


source,’vehicle is then tracked acoustically by hydrophones 
arrayed in a fixed resition on the oceén tcttom. The basic 


funeticn is that of reasuringe the transit time of the sound 


wave fre: the Source to different hydrorhones. Thase trarsit 
times eratle deterrineticon of distances ¢ 1 directZons of 


arrival which in turn are trianrulated to determine positicn, 
a series cf which defines the vehicle's track. 

As in mest enrineerins protlers, the requirements pleced 
on an accustic source used for tracking ere o*ten conflicting 
in na‘ure, with the end result being: en ercineering comtro- 
nise., Ver example, in order to achieve we bread pattern 
from a single elerent, the clercnt's dimensions must be 
snail relative to the wavelenrthof the sound transmitted, 
However, to cLtain the desired ccund pressure level, the 
¢lercnt miy have to be driven so hard that ecavitaticen (an 


undesirttle effect) results. Por many practical rea 
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therefore, it is essential that the directional character- 


istics of a proposea transcucer design be predicted prior 


to construction rather than measured afterwards. That is, a 


clear understanding of the relations tetween directivity and 


desirn constants is essential to properly design a trans- 


ducer for a spccific purpose. 


Study of the geometry involved in tracking an underwater 
venicle ‘ith ocean-floor mounted hydrophenes (assuming that 
‘he hydrophones aré not located at excessive depths) reveals 
that the slant range from the Vehicle to the receivers is 
normally several orders of marnitude greater than the dis- 
tance from the vehicle to the ocean floor, It is clear, 
then, chat most of the acoustic energy from a source mounted 
on the undertody of a vehicle shculd be directed obliquely, 
with a relative minirum being transmitted directly down at 
the ocean ficcr. Consideraticn of a large slant ranse 
situation dictates that the directionality pattern extend 
almest to the horizontal, however, not to the deuree that 
wlll result in surface reflectitns, 

A feeling fer the difficulty involved in practically 
achieving the radiation rattern previously described can 


he obtained ty studying the P.M. Morse [Ref. 2] solution 


‘for 6 radtelly vibrating strip which, extends-andefintialy 
foe “ gee oe Cae 


a roy . 
eee i ti FE 


along a cylinder, For the physical diviensians Gf the ’s 


B see 
bafiwe (vorpedo) and the frequency} (/Sk"2) Involved in‘ our 


case of interest, the Morse solusion ‘indicates the acoustic 
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will ensonify an area directly beneath the venicle. That 


is, the combined effect of the.cylindrical baffle and rela- 


tively hirh frequency focuses the acoustic radiation in a 


highly directional beam normal to the surface of the radic- 


tor. Obviously then, it can be anticipated that to achieve 


the obliquely oriented pattern desired, some method of 
countering this Sccusing must be incorporated in any pro- 
posed transducer desirn. 

The intent of this report then is to describe procedures 
for compu tine one characteristic of a sound source used for 


underwatér acoustic tracking - its acoustic raciation pattern. 


The an-rcoach will te mathematical in eecpe resulting in the 


develorrent of a computer rrorra:r. which will caleulate and 

plot .ne radieticn pettern of a flush mounted radiator located 
on the wall of a cylindrical shape, such as that of a torpedo. 
Tne stucy will enecr:iss three €2) proposed design confifrura- 


tions, céserited in the following section. 
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III, THEORY 


A. BACKGROUND 

The subject. of sound radiation from vibrating objects 
dean old one ana the Journal literature is replete with 
papers on the theory. it appears that the particular problem. 
at hand, radiation frem a finite element on a rigid cylin- . 
drical baffle has not received much attention. Morse and 
Ingard [Fef. 3] treat the problem of radiation from line 
sources on a cylinder. Laird and Coher [Ref. 1] developed 
a sOlution for the far field radiation oo a rectangular 
patch on a rigid cylinder, which ccmes closest to represen- 
ting, the present problem. Fer this reason, it forms the 
basis for the model of rectangular scurces and the b> trinine : 
point for the model of a cireular source, For the convenience 
of the reader, scme of the results of Laird and Cohen are 
summarized below. 

In their study, Laird and Cohen extended the theory of 
Morse (Ref. 2] to the case where the source 18 any separable 
function of the azimuthal and axial dimensions. In this 
thesis, their results will be applied to derive the equa- 
tions for the "Patch" and "Serment" sources. In addition, 
an extension of their work to the peneral case of a non~ 
separable source is invluded to describe the pattern of the 


"Disk" configuration. 


Ld 


The general approach developed by Laird and Cohen uses 
the cylindrical coordinate system shown in Figure 1. The 
method assumes that the source, which is mounted on a rigid 
cylinder, vibrates radially in such a manner that its ve- 
locity distribution may te represented as a separable func- 
tion of the azimuth and axial dimensions. Thet is, the 
velocity distribution has the same functional form in the 


azimuthal direction independent of the axial dimension and 


vice versa. This is shown diagrammatically in Figure 2. 
The boundary condition at the surface of the cylinder/source 


: can then be given by the following expression: 


o 


oF: ik g 
-dut , © eu : Are 
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: | (1) 
where the Fourler cosine series represents the azimuthal 
dependence and the Fourier integral represents the axial 
depéndence. The eobine series was arbitrarily chosen for 
mathematical convenience, 

Having established the boundary condititn, the general 
aspreseton for a combination of outgoing eylindrical waves 


of even azimuthal dependence, given by 


tk 2) 
Zz , 
ak, 


(2) 


: = -Lot . ss : ° : (1) 
pr(r,¢,z) = e = coe x Atk JH, (Kp de 
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is matched to the boundary condition at r = a (the surface 


of the cylinder). Introduction of the far field approxi- 


mation for the Hanke? function, conversion to spherical 


coordinates, and the solution of a Fourier integral by the 
method of staticnary phase result in the following general 
expression in spherical coordinates at a field point des- 


cribing the radiation from a source on a Pe. eylinder; 


( ; ei (kKR-ut ) F(k,, o ae 
p(R,6,o) = 2pcu, ———— x Chie 20smo 
= sin ov H'S2) (pasing) 


(3) 


To considey a specific source using this method, one 


needs only to specify its locaticn on the cylinder, physical 


dimensions and velceity distrituticn. Knowing the atove, 


the: Fouries cceftictonts, band descritinn the azimuthal deren 


dence. and the Pourler trensfort., CK)» Gescribing the axte 


deper.celice, Can 66 centeuzated., cnbetsbuticn of the -tourder 


eceffieilents and the Yourter transfom: into Fquation (3) 


results in an exrressicn cescribin;; the raclation for the 


particular scurce ecnsidered. The frequency dependence ‘ts 
incorf.orated through the wave nunter, "k", 

The derivation for the "Patch" configuration parallels 
the development for the case of the uniform reetengular 
source calculated "y Laird and Cohen and Lllustrates use 


of Equation (3) for determining the radiation from a specific 


separable ‘source. 
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B. PATCH CCNFIGURATION 

The "Fatch" array (see Figure 3) is composed of nine 
equal-dimensioned elements (each of angular width 2a and 
height 224) with the center element 180° out of phase with 
respect to the cther @lements. 

Consideration of Figure 4 shows that the source motions 
of the patch cannot be described by separable, independent 
functions of the cylindrical coordinates ¢ and Z. 

To avoid the complexities involved with a non-separablé 
source, the entire array will be viewed as a simple rectangu- 
lar cource with uniform distribution (see Figure 5). 

Likewise, the center element will te considercd as a 
Silipiey Unsfornly cxestce Pectancular source with c*rensiocns 
2a by eZee 

By subtracting twice the pattern function of the center 
elenent from that of the entire nine-element array, the 
desired recult is achieved. This assures the validity cf 
the linear superrosition prine’ple. 

The coefficients associated with the functions deseribing 
the source metion in azimuth are calculated using standare 


Fourier Series relationships, This results in: 
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For the ceriter element of dimenstons (20,22) 


m = Does saake 


For the entire arrey of dirensions (€a,€Z,) 
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For the center elerent of dinensions (24,2Z _) 


20: 
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For the entire array of dimensions (6a,6Z,) 
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Substitution into Equation (3) yields the following 


expression for the far-field pressrre due to the radiation: 


sin3k, Z., 
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pressure in fhe herisontel pliune for a sinrle element of 
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dimensions 2a by 22, 


pa 2pcU 2, ei (kKR-wt ) 
P(R,5,0) = 7 = — ——— x 


’ 
HS 


es esinro/male 
m=1 ea (ka) 


+ cosmé) (5) 

In order to sum the contribytions fror. all other elements 
which ere disycced at uniform anrles around the cylinder, 
it is necessary to transform the anrle (¢) in Equaticn (6), 
s0 that it reprecents ya the relative angle from the itr 
Clement to the fleld rcint, whee the fié14 point hee 
»f) measured fro*, the center of 


2 et 8 + te 
Srherica, seeriicctes (h, 
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qnene A, is the retative’ amp Lttude of the ith. elément’ and : ile 


vy is its relative phase and 9, = [ (4-193 + 8B). 9 


In a Simetee manner,. the, following result is eptained 


for pressure as a function of the angle in a plane con- 


G 


taining the axis of the cylinder: 
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D. CIRCULAR PISTON. SOURCE (DISK) gon aoe 
The "Disk" array is composed of a central cures 
piston source whieh is. Surrounded by 2 concentric annuler 


piston, with a 80 degree phase difference betweer, theix, \ 
{ My 
‘ 5 af of 


motions. | 


Since the motions of a circular piston aoe the. 


‘ 


side of a cylinder cannot be described as separable functions 


yu 


of the cylinder coordinates,.an extensicn of the méthods 
eseribed above must be made. Figure 7 shows (assuming a 
uniform velocity distribution) that although the functional 
form of the velocity aistribution remains the ‘game across 
the disk, the source is nonu-separable in that the axial 


limits are dependent on the azimuth dimension and vice versa. 
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Applying the boundary conditions for the circular 
piston using Equations (1) and (2) and the basic relation 


between particle velocity and pressure 
y 2-2 (8) 


results in: 


+% ik _Z 


Uniree * ~= en tut E Pada a k An (k, yur 2) (kale 2 ak 
(9) 


Equating: Equations (1) end (9) permits solution for A.(K5) 
in terms of the Fourier coefficients, an In the case of 
the separabie source, the Touricr coefficients, a,» are 
functions of the azimuth dimension of the source, ¢. Howe 
ever, in the case of the non-separable source, they wie aiso 
fiinctions of the axial dimension (2) and as such are expressible 
themselves as Yourier integrals. 

To COnpere the general Solution; having solved for 

Lpwl 2 Fk, bY 

A, (k 2)? (in the separable sade, A ks ie aTSLE TaN we 
Substitute into EG ation (2), the general Gallation . or 
outfoing cylindrical waves of even azimuth dependence 


which ylelds the following: 


p(r,o,z. Xr cosme raat 
: m=0 -0 (x, sa) 
, 4k Z 
C1)3. oe 
Cy Ckre * ale,) (10) 
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Introducing the far-field approximation for the Hankel 


function 
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#)) (xr) = oa 28 ge (11) 


results in: 
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Converting to spherical coordinates, results in: 
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Evaluating the Fourler Interral by the method cf stationary 
phase gives as a final expressicn for the radiation: 
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Returning to the non-separable disk, the boundary 


condition Equation (1), 1s again expressed. as: 


= 1 a7 lwt 2 oy 
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o ik,z 
f(z) = - F(k, Je dk, « 


Since each term of the sum over "m" 1s multiplied by 


f(z), we may inelude f(z) in the sum. ‘Thus, 
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Now, let By 62) = a (2) f(z), 


Define 
Gtk) = Pe (2d) just ag F(K,) = V{f(z)} 


Then Equation (16) becomes: 
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Let a, * 1 in Equation (1). Then Equation (17) is identical 
to Equation (1), which is from Laird and Cohén's develormen 


except that our GOK) depends on "m" while the F(k,) in 
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Equation (3) does not. However, th*s still permits using 


Laird and Cohen's results; namely Equation (13) still applies 


with a. in (13) set equal to one (1) ane G(k,.) replacing 
F(k,) ry 
The parallel to Equation (14) then is: © 


| ~ iis 

, et (kKR-ut ) eo G_(k, Je a 

p(R,6,6) = 2peU os) Ea zr 2 6b ES cosmo 1 
es o R sin 6 m=OH) 1 (kasinég) 


where we have set a, equal to one (1) in Equation (14) ane 
have replaced F(k,) by G(k,), but have included G.(k,) 
within the summation over "m" because Gk), unlike F(k,), 
depends on "nm", 

It remains, then, to evaluate 0(k,,) for the case at 
hand. The followine was aceorplished tr. collaboration with 
Steven R. Cohen [Pet. 07, 

Noting that the equation of the disk on the cylinder 
is given by. 2° + a? 6e zg where es atiy we Solve for 


"6" in terms of "z" 


Using, this result, we can express the Fourier coefficients 


Girectly in terns of (2), the axiel dimension: 


a, * 


9 
u 


2g 


eS ein ame ee 


Since f(z) = 1 (a constant) 
GK.) = Fla (z)f(2)} 


1 
= Pla(z)} = PES sin(B(z 292°)? 3} 


mn 
mf 2.2 mf, 2 .2 

a +o ta AG = a te¥0 a 

= - f Se ee e dz (20) 


To evaluate this intepral, we first let k, = w, 2 = t and 


z 
az = dt. Therefore 


nfo ee a? 
+e gm mgs al wife fo 
Gly.) = os f te No ~e NO ~ yettot oy 


o™ 
(21) 
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To evaluate the ebcve Tourier transorms we wlll use the 


follovine standard Leplece transform relationships: 


atVs*-a° 
(1) blot oen’)} = por tt > b 
(@) BY) ea naarE 


¢ 
Ss ta 


(2) L{c(t,a)} = Cs,n) 


(3) SL{C(t,a)}da = sC(s,a) da 


Applying relationship (3) to (1) 


153 (Ne 0" da} = s& da (23) 
s“ta° 
ate a s"4a" 
ab ‘° 
Letting t = in, wrere "c" is chosen vice "a" to avoid 


confusion, 1.e:, ¢ = a, ane substituting into (23) yields 


. ee : 1 ABs “42° 
L{fI (a toe )da} = - Sart © ‘ (2u) 
e a(-) 


RO ete Mette tha. Tetoesnt. Jt CEyt pe Ga 4 
2ur 2.m* 
Les us Weotes . Therefore du =Vt tos da, 
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Ce eee ee ee 
2 ee 
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ee } 
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Therefore Equation (24) beccemes 


(27) 


whieh is of the form Ff 


Unine the relaticne!.é: rér(t)) = anr(-w) with Equrtt 


(27) gives the result: 


Bie aes 
ce ~ 
F{- nin e . } 
! é 
= T fo vat eee) + J ( ) + feat 
"y ¢ ») 
a‘ . 
peice (2%) 
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Ff{e } 
me 
4 4 
oA a Yuet ) + a oFe tag 
2m ee 
a” (29) 
fon 
which can be directly used to evaluate 
th 
2 a o-” 
Pf mri e } 
where m= m, ¢ = a, and @ #® Zo° Therefore 
im 2 2 
Yo: Bie ee 
Fea. 8 } 
ri: 2 
1 © lin 4 \/ aes 
Sr Oa) fo, (2Q aoe) + J,( d+ oe] 
Bast % 
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soca via ve Ye a ay ) + Jenk ) + eed 
2 ooo an ao " 
\/ Oe 
aw se (30) 


im 
A similar tres¢rent lestine b = 7% permits the folleving: 
te ee “ 
eae an eg 
r 56 
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Since 


= Ff sin( tafe °- 
G (kK) F{-— sin( z Zz 


im 2 “=Z Amy, fi2 


Fit e8 9  } + Pre 2 © } 
G_(k_) = ef: ae eee ee: 2g) + tak. 24 ) 
m2 a “1 "ah &z ae 3h Mg Ne me 

Vk oe 
Zz 
a 
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z 2,m \/. 2.r' 
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Substitutine G(r.) into the reneral expressicn for the 
oe oa 
radtation, Equation (1°), results in the follovinr equation 


describing the racéiétion for a "disk" source en a eylinérical 


baffle: 
1G6pet,2° el(Reut) © @ 
p(t 6 aa 
a7? a hsané m=0 120 
aime 
Sa Gh FSET ETTS cosmé (33) 
He (kasinod) 
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The pattern for the "Disk" destin is achieved by sub- 
tractine twice the pattern of the inner circular piston 
source from that of a ciréular pistor source having a 


diameter equal to that of the outer concentric annular piston. 


aa aa ae ig: | 


III. COMPUTER MODELS 


A. GENERAL DESCRIFTION 

Based on the closed-form equations derived in the pre- 
vious section, a computer prorram was written to permit 
rapid calculation and plotting of the directivity patterns 
for the desirns of interest. A complete listing of the 
prorram developed is included as Appendix A. A brief des- 


eription cf each of the mafor subroutines which corprise 


Throurh the use of the prerrar, the parareters which 
effect the ractation rettern of a particular source can te 
yavieé ara thote ctfcet on the radisticn sattern noted. 

In thie runrer, a@ desicn for a rartieuler conflruraticn. can 
be achioved, 

che pataereters sPfeo tine the ypadictien pettern ere 
charred theceurh the use of “data inruts"” to the prorrar. 
The date inputs yormit the fcllowing: (1) ¢hanging the 
physteal Gimensicns of each of the sources, (f) changing 
the frequency at “hich the source is cperated, (3) changine 
the radius of the cylindrical baffle on which the source 
is mounted, and ¢"} 3electinr the plane in whieh the direc~ _ 
tivity pattern is @esired. In addition, for the use of 
the "Serment" desipn, the propram inputs permit specifying 


the amplitude and phase shadinr of each of the individual 
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elements of the array. However, in this case only the 
resulting bear rattern in the horizontal plane can be 
plotted. 

Another input which was added to increase the flexibility 
of operation of the program is the “summations limit". As 
can be noted from the equations derived in the previous sec- 
tion, all of the radiation pattern equations include an 
infinite series summation. Althcugh this limit remains 
an input variable tc the prorram, judicious use of this 
input 18 recorrended 1 correct results are to be obtained. 
This ratter will be discussed further in the "Results" 
section cf this thesis. 

Proprar cutruts are both tabular and rraphicali in 
character erd can be summarized as follows: (1) The cor- 
puted values of the rest and imarinary parts of the cornlex 
pressure fcr cre (1) derrec tearing tnerements expressed 
in dynes per square :éter, (2) The rarnitude of the cer- 
plex pressure for one (1) derree bearing increments expressed 
in dynes per square meter, (2) The marnitude of the complex 
pressure for conc (1) derree bearing increments expressed 
in decibels refererced to one microtar, (4) A polar plot 
of the rmopnitude cf the cerrlex pressure 
norralized to the larrest value cf the complex pressure, 

(5) A polar plet of the somd pressure level in decibels, 
normilived to the lével of the larrest lobe, using a scale 


lenrth of 50 Ak, 


3( 
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B. DESCRIPTION OF COMPUTER PROGRAM OPERATING PROCEDURES 
A sample input data deck is included as Appendix C to 
illustrate these written instructions. 
With the exception of the amplitude and phase shading 
portion of the "Segment" design, the input data decks re- 


quired to ottain directionality patterns for each of the 


aan prepmcnonges oe 


three (3) configurations are identical. The first data 
card indicates the configuraticn for which e beam pattern 
is desired. This is accomplished by typing the name of the 
configuration (PATCE, DISK, OR SEGENENT) as shown on the 
‘irst data card commencing in column number one (1). The 


second date cerd indicates the number of summations desired 


ticular desirr for which the pattern is being calculated. 
The nurter of surrations desired is typed in columns one (1) 
and twe '2), [FCRMAT I27, oneshalf the height (25) of a 
single eiement measured in reters is entered in columns 
11-20 [FOPMAT F10.5) (in the ease of the DISK the outer 
radius measumeé in meters), one-half the angular width of 
a single element (a) measured in radians is entered in 
columns 21-30 [FORMAT F10.5] (in the case of the DISK - the 
inner radius reasured in meters), the radius of the cylinder 
| measured inmaters in columns 31-40, ard frequency in KHz (i.e. 
for 75 KHz type "75.") in columns. 41-59 [FORMAT F10.5]. 


The third and final data card indicates the plane in which 
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a pattern is desired and the angle which is to be held 
constant for the calculation. By reference to Figure 1, 
one can see that by setting @ equal to a constant and sum 
ming 6, a pattern in a plane containing the axis of the 
cylinder is achieved. Likewise, by setting 6 equal to a 
constant anc summing ¢, a pattern for a plane orthogonal 

to the axis of the cylinder is obtained. In this manner, 
various planes for a particular source and three (3) dimen- 
Sional aspects of the pattern can be perceived and investi- 
gated. Indication of the angle which it is desired to vary 
is indicated by typing a "1" (Theta varying) or a "2" (Phi 


varying) in columns 1-2 [FORMAT I?]. The angle held constant 


is indicated by typing the anrie (in derrees) in columns 
11-20 [FORMAT F10.5]. 

To obtain a pattern jn the horizontal plane for ampli- 
tude and phase shading for the "Segment" design, enter a 
"1" or "2" in column “2" [(PTORVAT I2] of the fourth card of 
the Seprent Data Deck. On the subsequent card (the fifth 
care of the data deck) enter the total number of elements 
to be shaded using columns 1-3 [FORMAT 13] (4.e. if 24 
elements. are to be shaded, tyre "24" in columns 2 and 3). 
Followinr this care, a separate card for each of the ele=+ 
rents requiring shading is entered. In columns 1-3 the 
elerent numter to be shaded is entered (FORMAT I3], in 


cOlumns 11-20 the arplitude of the shading is entered, 


(FORMAT F1LI.6] in columns 21+70 the desired phase shading i 


entered [(!ONVAT F10.5] as an anpyle, 
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‘IV. ‘RESULTS 


A. PATTERN FOR UNIFORMLY EXCITED RECTANGULAR PATCH 

Figure 8 shows the horizcntal directivity pattern calcu- 
lated for a Single element (a = 3.7°, ka = 14) using the 
"Segment" subroutine. The pattern was calculated and plotted 
in rectangular coordinates to permit a direct comparison 
to a similar pattern plotted by Laird and Cohen [Ref. 1] 
using the same input variables. Comparison of the two 
patterns shows they are identical. This result verifies 


that the basic coding of the program is correct. 


B, EFFECT OF VARYING THE. NUMBER OF SUMMATIONS 
As noted previously, judicious use of the eunming 
limit is required if accurate results are to be obtained. 
@ - & Runs for the same configuration were mace with N = 5, 10, 

A A 20, 25, 30, and 35. Little to no observable charge was 
noted in patterns for ‘the Syne with "HN" greater than twenty 
(20). However, fon runs with "Nf less, than twenty (20); the 
pattern varied considerably for each run, Experience with 
the program indicates that about 20-25 terms of the infinite 
‘serics (which eppéars in the pattern function. of all de- 
signs aonetdoned) must be calculated before the diverrence 
of the Hankel derivations in the denominator cause the solu- 

az tion to converge. This conéelusion concurs with that errived 


at by Laird and Cohen in their original development. 
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As indicated by the comments in the BESJ SUBROUTINE 
listing, a maximum of twenty summations should be calcu- 
lated if the value of the entering argument is less than 
fifteen (15). 

To meet both the previously mentioned requirements, 
it is recommendéa that the summing limit always be set 
equal to twenty (20) for any final computations of a pro- 
posed design. The ability te change the number of summa- 
tions has been included to permit a rapid rough first-cut. 
The execution time for the eccsputer solution varies from 
5-20 minutes cependent on the number of summations and the 


input parameters, 


Ce FREQUBNGY URE beaCh CF Tee SOLE: 

Figur®s 9, 10, I1 and 12 are included as examples of 
the manner in which the frequency dependence of the patterns 
can be Studied. All input variables were held constant 

\ 

and the frequency vas varied to obtain these patterns! The . 
important desirn capability provided by this feature is the 
ability to predict the effect of the frequency bandwidth — 


of a proposed design on the acovstié radiation pattern. 
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Ve. CONCLUSIONS AND PECOMMENDATIONS 


Computer models which permit the computation of the 
sound radiation pattern for three (3) different source 
configurations mounted on a rigid cylindrica? baffle have 
been develcped. 

The computer solutions agree favorably with rrevious 
mathematical results ottained by earlier investirators. 
Althourh erpirical dita are not availetle for confirmaticn 
of the patterns, 1¢ Js considered that rood apreement would 
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APPENDIX B 


MAIN 


Purpose: To control prorrar 

Method: The program inspects the first four (4) letters 
of the ecded tities (DISK, PATCH, SEVEMENT - Note the 
spelling of the "Cerment" desirn) to transfer control to 
the rreper subroutine. It also writes out the first line 
of output idersifvine the desirn confircuration being 
ealculate?, 


pas sia + 7 sy Re ee ; 
Calvec hy s> Noyes ins: eaea Care 
FLAS Re Re MEL ae 


eotre 
Cay e Sey Ss 
Sasa ceed ee. 
tay rerensirs 
ti ao 
bees 


oo | oo * o*. t. . © yee, tae . 7 
Porpcsey  Teoooeeteetate ng Nanrol derivative piven the order 
Lu 


VOSS: a SUP ee Set OMe Or ero wear ae nl 
COLOT ESE Oe Be aie es bao nee poeur ne poe relat cs, 4 
REC) © Ee Bay Gs. She Henkel Sunetion of 

(8) : ee (iby) : 
order! SE) ed ication OCP) eee ap yentatee te eorbiaing che 
WD gthegh A OR SR AR tetas. tay hs ie bee Ne 


2 oT Tree 
Called fv: TICE 
te oe eat : 
SALUT 
PAUSCH 


SEGEENT 


alis Te: RES 
BrSd 


Te 


DISK 

Purpose: To calculate the radiation pattern for the "Disk" 
design. 

Method: The "DO" loor ending with statement number nine 
(9) calculates the pattern for the "DISK" design. It con- 
tains within it, a second "DO" loop ["DO loop ending with 
statement number eirht (&8)] which calculates the infinite 
series ccrrosed of the odd Bessel Functions. The entire 
loop is executed two tires through the use of the indicator, 
"¥MLAG." After the seccnd executicn of the loop, twice the 
results of the first execution are subtracted from these of 


the secend execution. The statement, "IF(CAFLAG.EG.1) Gd 


TO V4" transfers the results of the calculations to the 


Sica emer Conk ete 
POA ST 
PRECH 
Piuvsnces Ce entoutete the Yosletien p atten for the “rete” 


design. 


PelVOns \ The pusi@-eruny ion’ denived for the. "Pateh desiin 


are aerded tn ¢he "Po" Jeon ending with statement number nine 
(9). The pattern is achleved ty exeenting that "DO" icc} 
twice. In the first execution of the loop, the pattern for 
the oppositely phase shaded eenter clement is calcilated. 


Subsequently, "XFLAG" (an indicator of how many times the 
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rT 


ate 


loop has been executed) is updated, and the dimensions of 
| the element changed to those of the outer (larger) element. 


After the second execution of the loop, two (2) times the 


[ first pattern is subtracted from the pattern obtained during 
{ 

E the second executicn of the loop. The "IF" statement, 

| IF(KFLAG.ES.1) GO TO 10" transfers the calculations to the 

| plotting subroutines, 

| Called Fy: MAIN 

t Calis Te: HANK 

POLPR? 

I 


SEG 

5 m : Ion ° 
Furnerns: Ts Galoulwt>e. fhe: redieticn. setter. for the 
pais LeeeaS 


eo ns ” 52 2 es 
BOOPUERL! est on, 


eo Ir er pe eo 


endine vith cteteras 4 nurker efy (6) varies cither anrle 4 


ne f a 7 r) A , + ers + ; 6 
rethors Chak SUIS re coy. Pe. Shouse 26 Consist “or “td 
ce 
t ae 4 bens 4 ry 4 i] say ba il 2 
PAPess Fi) Stator te iy revtous, cote Staten ert. MRELOOCr 27) 
I 
{ wy ff ’ 5 . Sepia yee Ye pa i oF een 
I His Pt. Oe ee 2s Serra ie: Parser: fOr a Sie he 
‘ * 2 sS 2 % Pa 4 bal ~ 
CLEMO igs hes Oo A Td SORES OMA’ ee at Me Se EMO (ORTON ANE 
} : : , 
3 , . ne * ae 
i Che peor Pes cee: Bee Ts Pees 2. Toe Oye. Te. 
t 7 s a . “ * 
eG eben peese ree ars Ae, oe eee Searle, ee eo are a aes Fe aie ae Gs 
h 
Ne - 4 "é@ te a tery ee 3 an ne 4 
Deperder* ety “he: Dee On Fee Ee he SOE 
i 
t 
Ly 
[ 
& 
k 


OP “atipie ¢ “pore ly ep fe Dre Sees, 
The anple held eonstant (second entry on the third card 
Of the inrut data deck) in consonance with the anrle te be 
varied (dictated Ly the "IPL" code) specifies the plane in 


which the pattern is determined. 


i 4 


In the second part of the subroutine (amplitude and 
phase shadinr), the "DO" loop ending with statement number 
sixteen (1€) is the coded version of the general expression 
for the amplitude and phase shading derived earlier for 
the horizontal plane. 

Called Py: MAIN 


Calls To: HANK 
POLPRT 


POLPFT 
Purrere: To control the plotting of the polar plot. 


Method: This subreutiné fs the main subroutine itn the reler 


Plot pecwace end is resrensitic for calling the varicus 


ro ie . € @¢Se . Yar ~ 
tones cf the pac“laore, 


7) 
Cc 
t oad 
es. 
UP 
= 


The seale facter, S, must be chanred aceerdinge to the 
printer chiracsterictics. The seale factcr in this subreutise 
is set. for -ten, Cy eharacters per inch: Sor the etseisse 


find eirhe *, shennceters per inch form-the oreinate: aria, 


After 2net sfc l eine Te ai, veces me &, the Jurus: date, 
Ys 2h S¢arees te cclerm se re sorrpelivive feetor. 2f this 
normsilivine fretcr 2s lees thin l.he@?, an error statement 
TecT yon en Ce te ene Te eLertad. 

In the LO LOG! eiding with statement &, each line of 
the polar plot is printed after a call is made to PRPLOT to 


establich the polar jirrid information. The variable, DIM, 


4s used as a scaling factor for the polar plot. The value 
of 1.0 will cause all of input data to be plotted, however, 
4f only the values less than one-half of the normalizing 
factor are of interest, then DIM can be set to .5. This 
will enlarge the center of the polar plot. 

DISK 


Calls To: PTPLOT 
SAPT 


PTPLOT 

Purpose: To establish the grid informaticn for the polar 
plot. 

Methoc: In the DO LOOP ending at steternent 1 the aiprie 
numeric characters arc transferred to ISY! in order to pass 
via COMMS to other cubroutines. In the etatenents follevine 
statement 2, the eauetiens for the rictted concentric cireles 
are esteblished. Relcu ststerent 7 the crid marks on the 
090-270 axis are iIncerted. 


Called Py: POLPET 


al T avons, 

Calls Ta: LIMECL 
yrry bad) 

as hidas 
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<n rere TE 


ll Eo 
: 
4 


NUMB 

Purpose: To place derree numbers on the polar plot. 
Method: The current line which is being printed is passed 
to the subroutine in the callinr argument. If this line 


ecortains derree nurters, these numbers are placed in the 


correct rosition ty the IF statements. 


Cal.8: Sex HONE 

eee Sats 

Furpese; To insert prid characters cn the polar plot. 
“Nothods The pertod charteter (ITSi(2)) is inserted in the 
prever ree tice. in the statetente ateve stecerent ty. “in 
the statements after statement 4, the rrid nurbers latels 


Callet bys POLPECT 


ade ste 
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